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The invention is described in the following statement 
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SULPHUR DIOXIDE DETECTION METHOD 



Field of l:lie Invention 

5 The invention relates to a sulphur dioxide 

detection method and apparatus. 

Background to the Invention 

10 Volcanic ash and sulphur dioxide clouds 

constitute a serious hassard to aircraft even after the 
clouds have rooved from the site of a volcanic eruption. 
Apart from containing ash particles, the clouds include 
gases such as SO2 which after a few days oxidises and 

15 hydrolises to form sulphuric acid droplets, either as an 

aah-acid mixt\ire or as a coating over ash particles. Both 
the ash particles and the sulphuric acid droplets of 
volcanic ash clouds are capable of causing significant 
damage to and possible loss of an aircraft which 

20 encounters am ash cloud. 

A nuBiber of aircraft encounters with volcanic ash 
clouds or sulphur dioxide clouds have been recorded in the 
past where significant damage has occurred. It will be 

25 appreciated that the sulphur dioxide may be found in areas 
separate from the volcanic ash. In the year 2000, a 
National Aeronautics and Space Administration (NASA) DC- 8 
Airborne Sciences research a irpl ane flew through what was 
described as a diffuse volcanic ash cloud from the mount 

30 HEKLA Volcano when flying from Edwards, California to 

Kiruna, Sweden. The ash cloud was not visible to flight 
crew, however, the science research airplane carried 
sensitive research equipment which was capable of 
detecting the sulphur dioxide. In-flight checks and 

35 post-flight visual inspections revealed no damage to the 
airplane. However, detailed examination of the engines 
revealed damage to some of the turbine cooling passages. 
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Furthermore, high levels of sulphur were found in the oil. 

It seems likely that this ash cloud actually was 
predominantly a sulphur dioxide aloud. Even if it was 
5 not, it raises the possibility that an aircraft can fly 
through sulphur dioxide without passing through ash. The 
post encounter treatment of the engine in the case of 
sulphur dioxide encounter would be different to and 
considerably cheaper than the equivalent treatment 
10 required of an engine during an ash encounter. 

Accordingly, it would be desirable to provide a 
sulphur dioxide detection technique. 

15 Summary of the Invention 

The present invention relates to a method of 
detecting sulphur dioxide clouds comprising: 

measuring Infrared radiation at a viewing 
elevation at or above the horizon and at a Key wavelength 
at which there le an SO. feature and in the vicinity of 
which there is a region where the amount of infrared 
radiation from water vapour in the atmosphere varies in 
accordance with a predetermined relationship; 

measuring radiation at two or more subsidiary 
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wavelengths in said region; 

determining the amount of radiation from water 
vaBOur_at the key wavelength from the measured radiation 
at the subsidiary wavelengths using the predetermined 

30 relationship; 

determining whether a sulphur dioxide cloud is 
present from the measured infrared radiation at the key 
wavelength and the determined amount of radiation from 
water vapour. 

35 

Preferably, said subsidiary vavelengtbs are 
located either side of said key wavelength. 
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The inventor has determined that the key 
wavelength should be selected from the group of 4jjm, 
V.Sjim, and 8.6faBand that 7.3Mm is the preferred key 
5 wavelength. 

Where the key wavelength is T .2]m, it is 
preferred that subsidiary wavelengths at iO.Sjun are used. 
The inventor has established that for the region of these 

10 wavelengths the predetermined relationship is that 

radiation from water vapour varies in a siibstantially 
linear manner. Accordingly, the radiation from water 
vapour at the key wavelength can be interpolated from the 
radiation at the subsidiary wavelengths on the basis of 

15 this predetermined relationship. The inventor has also 

established that there is substantially less SO, absorption 
at this wavelength. 

The method may also involve cos^ensating for 
20 background SO2 in the atmosphere. 

The invention also provides a detection apparatus 
for detecting a sulphur dioxide cloud comprising: 

measurement means for measuring infrared 
25 radiation at a viewing elevation at just below, or above 
the horizon a key wavelength at which there is an SO, 
feature and in the vicinity of which there is a region 
where. the amount of infrared radiation from water vapour 
in the atmosphere varies in accordance with a 
30 predetermined relationship and for measuring infrared 
radiation at two or more subsidiary wavelengths in the 
region; and 

processing means for determining the amount of 
radiation from water vapour at the key wavelength from the 
35 measured radiation at the subsidiary wavelengths using the 
predetermined relationship and determining whether a 
sulphur dioxide cloud is present from the measured 
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infrared radiation at the key wavelength and the 
determined amount of radiation from water vapour; and 
output means for generating an output signal 
indicative of the presence of a sulphur dioxide cloud when 
a sulphur dioxide cloud is present. 

The inventor has also determined that the method 
and apparatus of the present invention can be used to 
detect sulphur dioxide clouds from the ground or from an 
aircraft. 

Brief Description of the Drawings 

Figure 1 illustrates the S02 absorption feature 
in the region 1200cm-» to ISOOcm*' and the preferred 
measurement wavelengths of the invention; 

Figure 2 is a schematic diagram of a SO2 
detection apparatus of the preferred embodiment; 

Figure 3 illustrates two modes of operation of 

the apparatus; 

Figure 4 is a schematic diagram of apparatus to 

be used from an aircraft; 

Figures 5a - 5c represent normal climatic 

conditions; • 

Figures 6a and 6b represent variations on normal 

conditions to allow testing of the invention; 

Figure 7 represents variations in SO2 for 

testtjRgi .-^ on' 

Figure 8 shows variation in temperature with SO2 

concentration; and , 

Figure 9 shows temperature plotted as a function 

of absorber amount. 

DescrlptiPn of the Preferred gm>.«,<<inftnt 

Herein, the term «key wavelength" Is used to 
refer to a wavelength at which there Is an appropriate SO, 
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Persons skilled in the art will appreciate that a 
^wavelength" in the context of this specification does not 
iinply a single wavelength but rather encoa^passes a band of 
radiation. Typically the width of the band will depend on 
the filter used to observe/measure light at the wavelength 
of interest. 

The term "subsidiary wavelength" is used to refer 
to a wavelength in a region in the vicinity of the key 
wavelength where a relationship can be established between 
radiation from water vapour at two or more subsidiary 
wavelengths and radiation from water vapbur at the key 
15 wavelength. 

The preferred enibodiment provides a method and 
apparatus that allows identification of sulphur dioxide 
clouds in the free atmosphere. The apparatus of the 
preferred embodiment uses an infrared detector, 
interference filters and focussing optics. The filters 
divide radiation within the band between 6.8 and 8.1 H» 
into three narrow bands. The central bands corresponds to 
a strong SO, absorption feature caused by the antx- 
symmetric stretch of the SO3 molecule. The other bands are 
above and below this feature. The central band B„, is 
sensitive to SO, concentrations. The lower band, B. «xd 
Higher, band B..,. are used to account f or . the_ef f^cts of 
water vapour on the absorption in band Bo- 

Accordingly, Be is the key wavelength and B, and 
Bh are the subsidiary wavelengths in the preferred 
0xiibo<iixnent * 

,5 Figure 1 Illustrates the absorption feature due 

to SO3 for the infrared region extending from 1200 cm 
(8.33 H»> to 1500 cm-^ (6.67 Mm). The ordinate in this 
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plot ia line strength and the abscissa is wavenumber (cm"*; 
wavelength in fun » 10 , 000/waven\iinbor in cm*,^) . Also, shown 
are three idealised filter response fiinctions which 
isolate radiation within the three narrow regions 
corresponding to: B* (7.633-8.065 jm) Be (7.143-7.57 }m.) 
and Bi (6.897-7.042 urn) . 

The response functions are noiroalised to unity 
and scaled appropriately for plotting. The central 
wavenuniber for the SO2 absorption is 1363 cnT* and the band 
extends from about 1320 cm'* to about 1390 cm**. A filter 
covering this region responds to all the radiation from 
this band; whether the SO2 feature be due to absorption or 
emission. In the case of a detection apparatus viewing a 
cold background, i.e. viewing from the ground to space or 
from an aircraft towards the horizon, there would be more 
radiation in this band in the presence of the SOa cloud 
than if it were absent. 

In practice, water vapour and clouds also absorb 
and emit radiation in the region 7-8 ^m. The inventor has 
realised that the two bands positioned either side of the 
central band are used to eliminate the effects of water 
vapour . 

Water vapour absorbs and emits radiation 
throughout the region 7-8 juft. The amount of radiation 
. . ..atesojebed or emitted .depends .on_ t^^^ 

and on its location in the atmospheric column. Water 
vapour near the boundary of the earth's surface is 
generally warm and abundant. Water vapour near the 
tropopause (i.e. at jet aircraft cruising altitudes) is 
cold and sparse. The central band of the SO3 detector 
responds to radiation due to both SO, and water vapour. 
The lower and higher bands of the detector however, are 
only sensitive to water vapour. The inventor has 
determined that the radiation from water vapour in the 
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region surrounding Be behaves in a sufficiently linear 
manner to enable it to eliminate the effects of water 
vapour on the central band. The Planck blacJcbody 
radiation from Bi and Bi, are linearly interpolated to 
estimate the radiation detected in Bo due to water vapour 
only. This radiation amount is subtracted from the 
radiation actually measured by Be. The residual amount is 
due to SOa* 

A schematic of the detection apparatus is shown 
for illustrative purposes in Figure 2. The detection 
apparatus 6 consists of four major components: 

• Fore-optics 1 that focus a beam of incoming infrared 
radiation onto a detector. 

• A filter wheel 2 consisting of at least three narrow 
band interference filters that isolate radiation into 
the bands I Bi, Bo and Bj,. 

• An infrared detector array 3 sensitive to radiation in 
the 7-8 |4m region. 

• Processing means 4 for processing the detector signal 
to determine whether SO, and hence a sulphur dioxide 
cloud is present. 

Figure 3 is a schematic illustrating two modes of 
operation of a detection apparatus 6 that senses infrared 
radiation in order to detect SOa clouds. One mode assumes 
that the detection apparatus 6 is on board an aircraft 7 
an d vi ews the SQa cloud !BL)Mad_at a mall to the 

horizontal. The second mode assumes that the detection 
apparatus 6 is baaed on the ground and views the cloud at 
a large angle to the horizontal <e.g. zenith viewing). 

The detection apparatus of the preferred 
enibodiment may be operated from the ground viewing the sky 
above or^from an aircraft viewing forwards- The principal 
mode of operation is anticipated to be from an aircraft 
with the instrument having an unobstructed view of the 
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atmosphere ahead of the aircraft as the method works best 
when water vapour concentration is less than Ig cm . 
Ideally the view should be horizontal or a few degrees (3- 
5») above the horizon, so that the background radiation is 
cold. Typically, aircraft fly with their nose at an angle 
of about 3 degrees to horizontal. However, the processor 
4 can be configured to account for changes in viewing 
zenith angle, making the technique insensitive to the 
viewing direction. For the case of a detection apparatus 
6 viewing ahead of an aircraft at a zenith angle of Z 
degrees, the detection apparatus 6 provide.s three signals 
to the processor 4. A synthetic signal corresponding to 
the amount of radiation from water vapour is determined 
through linear interpolation of the signals from Bi and B^. 
This signal labeled B, is compared to the signal from Be- 

The processor 4 then cons>utes the SO3 amount at 
the key wavelength Be using 4 and the original signal B„. 
The processor 4 uses pre-defined look-up tables that 
account for standard atmospheric conditions (tropical, 
mid-latitude, and polar) and the viewing geometry. The 
detector array 3 provides an image of the SOa amount with a 
spatial resolution that depends on the exact number of 
detector elements in the array (320x240 is recommended) 
and the distance to the SO, cloud. Distance information 
cannot be supplied by the detection apparatus 6, however, 
the SOa anomaly will be detected at distances of up to 
several l00..3ms._4SP.ending on the .cruising altitude and 
clarity of the atmosphere ahead. The detection apparatus 
6 produces an output 6 either in the form of an amount of 
SOa or an alarm signal indicating the presence of SO,. The 
alarm signal may cause an audible or visual alarm in an 
aircraft. . 

Figure 4 illustrates how the apparatus works in 
the case of being mounted in an aircraft. 
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In addition to signals from the detector 3 the 
processor 4 also receives aircraft altitude Inforxnation 8 
from the aircraft and standard atmosphere information 9 
from a memory associated with the processor. 

Examples 

A sophisticated radiative transfer model-MODTRAN 
(Berk, et al . , 1989) is used to model the response 
expected from a single -element detector viewing aroalistic 
atmosphere. The viewing geometry is varied in the 
simulations to account for viewing from below the SO2 
cloud, viewing from above, and viewing at a small angle 
along a nearly horizontal path. The amount of SOa is 
varied, as is the main other gaseous absorber in the 
region-water vapour. We refer to the amount of SO3 as the 
cloud thickness. 

1. Model Atmosphere 

Vertical profiles of the model atmosphere used in 
the simulations are shown in Figure 5 and variations used 
to test the present invention are shown in Figure 6 and 
Figure 7. 

(a) Temperature 

The ten^erature profile is shown in Figure 6a. 
Varying the profile" has little effect on the retrieval and 
detection algorithm because the algorithm uses differences 
in teinperatures. No further simulations were performed on 
this parameter because of its insensitivity. 

(b) Watez' vapour 

water vapour was varied by increasing the amounts 
in the lowest layers from less than O.l cm of precipitable 
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water to more than 3 em. No effect was found on the 
detection or retrieval because the water vapour lies below 
the SO2 cloud. Water vapour was also increased in the 
layer that contained the SOj and this has a major effect. 
The perturbed water vapour profile is shown in Figure 6b. 

(a) Sulphur dioxide 

The vertical profile of the background SOa is 
taken from the US standard atmosphere. The profile 
corresponds to a well-mixed gas with a constant vertical 
concentration of lO"' ppmV (parts per million by volmne) . 
Perturbed profiles, with increasing SOj concentration, are 
shown in Figure 7. Ei^ht profiles are shown. The 
integrated amount of SO, in a vertical column for the 
profiles varies from 10 milli atm-cm to XOO xailli atm-cm.. 
Depending on the patblength travelled the total absorber 
amount can be much larger. Results for SO, absorber 
amounts of more than 1000 milli atm-cm. 

2. Viewing the SO. cloud al onff horizontal paths 

For the purpose of example, model simulations have been 
performed for the case of horizontal viewing from a 
platform (e.g. an aircraft) directly ahead and towards an 
SOa cloud. The viewing direction is assumed to be 
horizontal at the altitude of the platform (8 km, or 

-26,000 feet is assumed). «ie _?\«?»d thickness (as 

measured in the viewing direction) is varied from 10 km to 
500 km and the concentration within the cloud is varied 
from background levels to «0.1 ppmV. This range of 
concentration covers the smallest eruptions (that are 
likely to reach these heights, e.g. Hekla-style eruptions) 
to the largest observed this century (e.g. Pinatubo- style 
eruptions) . The results of these model simulations are 
summarised in two figures. Figure 8 shows the variation 
of the temperature anomaly (the temperature difference 
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between the synthetic signal and the measured signal as a 
function of cloud thickness) . 

The faxoily of curves 20-27 generated from the. 
modelling are lines of constant concentration for SO2 
concentration varying from 0.0136 ppmV 20 to 0.1083 ppmV 
27 . The points that lie on vertical lines correspond to 
lines of constant cloud thickness. As the cloud thickens 
the curves follpw the same trend with increasing anomalous 
signal until. the cloud starts to become opaque. At this 
point, which varies with SOa concentration, the temperature 
anomaly increases towards a limiting value (AT « -2 K) . 
Note that the opaque limit is reached either by increasing 
concentration or increasing cloud thickness, since both 
quantities increase optical depth and hence absorption. 
Beyond a thickness of 500 km, the cloud is essentially 
opaque and the radiative process changes from absorption 
to emission. 

Figure 9 provides an alternate way of 
widerstanding the physical processes involved in SO, 
detection. Here the texnperature anomaly is plotted as a 
function of absorber amount. The plot tells us that for a 
given anomaly, several values of absorber amount are 
possible, depending on the cloud thickness and 
concentration. Thus, it is not possible to uniquely 
quantify the absorber amount from the temperature anomaly 
without knpwlng either the con^cent^^ . 
thickness, in practice it is not necessary to know these 
quantities, as the purpose of the invention is to detect 
the presence of SO2 gas in the free atmosphere, rather than 
quantify the amount. The modelling does give an 
indication of the limits within which detectibn of SO, is 
possible. At the lower end, for cloud thicknesses of 10 
km or less, the SO2 concentration must be larger than »0.06 
ppmV. This corresponds to an absorber amount of «25 milli 
atm-cm. SO2 clouds that intercept air-routes (i.e. heights 
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>20,000 feet) will have horizontal dimeneions of 10' s of 
kilometres and absorber amounts well in excess of 25 milli 
a tin- cm would be expected* 
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